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I . INTRODUCTION 

Methods for studying the effect of variations in param- 
eters of a multivariable, multiloop control system have been 
developed by Mitrovic, Siljak, Thaler, et al (1) , (2) , (3) . 

Mitrovic's method consists of specifying as variables 
the two lowest order coefficients of the characteristic equa- 
tion. Using the characteristic equation as a mapping func- 
tion, constant zeta, omega, and sigma curves are transformed 
from the s-plane into the variable coefficient plane. This 
technique permits the adjustment of these variables so that 
the roots of the characteristic equation may be set at de- 
sired locations. 

In 1964, Siljac extended Mitrovic's Method into the 
Coefficient Plane Method (2). This technique permits the 
utilization of arbitrary pairs of coefficients as variables. 

The Parameter Plane Method is an extension of the above 
methods, which permits the adjustment of system parameters 
which are located in more than two coefficients of the char- 
acteristic equation (5), and (6). The Parameter Plane Method 
permits the analysis of complex systems which mav be either 
linear or nonlinear. However, the complexity of this tech- 
nique requires the computation of the curves by a computer. 

In 1965, Nutting developed a computer program (PARAM-A) 
which calculates and plots the Parameter Plane curves for con- 
stant zeta, natural frequency and sigma (7). In general, 
analysis of the resulting curves results in the location of 
all of the roots of a characteristic equation for a given 
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parameter (a, 3) pair. However in 1966, Bowie discovered sys- 
tems for which the existing Parameter Plane techniques did 
not locate all of the roots (8). Bowie designated these sys- 
tems as singular line systems. 

Section III describes the mathematical basis of singular 
lines and methods for proving their existence. A method is 
also developed for obtaining singular lines. Section IV con- 
tains examples of singular line systems. The application and 
analysis of singular line techniques to specific problem 
areas is contained in Section V. Section VI concludes this 
report with comments on existing singular line systems and 
recommendations for further investigation of singular line 
systems . 
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II. THE PARAMETER PLANE 



In this section the derivation of the parameter plane 
equations for the linear case will be derived (9) and gen- 
eral rules for mapping points from the s-plane into the 
parameter plane will be discussed (3). 



A. PARAMETER PLANE EQUATIONS 



Consider the characteristic equation: 



n 



F ( s ) = l a, s = 0 



( 2 - 1 ) 



k=0 



where : 

the coefficients a^ (k=0 , 1 , . . . ,n) are real 
functions of two system parameters desig- 
nated as a and B and s is the complex var- 
iable. Let s be written as: 

s = + j w n *^l - ? 2 (2-2) 

where w n is the undamped natural frequency and £ is the 

damping ratio. Letting 

0 = cos ^(-£) (2-3) 

and substituting into the equation for s gives: 

s = a) (cos9 + j sin0) = to e J . (2-4) 

n n 

Then s may be written as: 

s^ = u^e^® = u)^(cos k0 + jsin k9) . (2-5) 

n n 
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Defining : 



T k (-C) 


= cos k6 


= cos [k cos 1 (-?) ] 




U k ( -«) 


_ sin k0 


sin [k cos ■*"(-?) ] 


(2-6) 


sin0 


sin0 


And substituting 


into the 


equation for s results 


in : 


k 

s = 


to k [T, (-?) 
n k 


+ j/l-? 2 U k (-5) ] • 


(2-7) 


It may be shown 


that : 








v-f - 


(-l) k T k (?) 






u k (-c) = 


(-D k u k (?) 


(2-8) 


where T k (?) and 


U k (?) are 


Chebyshev functions of 


the first 


and second kinds 


respectively and are given by the recursion 


formulae : 








T k + i (5) - 2T k (t) 


+ T k-1 (?) = o 




u k + i (!) - 2u k (tl 


+ U k _! ( C ) = 0 


(2-9) 


where : 








T 


o ( ?) 1 


U Q (?) = 0 




T 


x (c) = C 


u x U) = -1 . 





Substituting equation (2-9) into equation (2-1) and equa- 
ting the real and imaginary parts independently to zero re- 
sults in: 



n 



I 

k=0 



a, co k (-l) E T. (?) 
k n k s 



0 



L Vn<- 1)k+ \ U) = 0 



( 2 - 10 ) 
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From equation (2-9) the interrelation between the Chebyshev 
functions is obtained: 



T k ( 0 = - U k _ 1 (0 . ( 2 - 11 ) 

Substituting into equation (2-10) : 

Jo ( ' 1 )ka k"n U k-l (?) = 0 

(2-12) 

I (-1)ka k a3 n U k (?) = 0 • 
k=0 k n K 



Assume for the linear case that the coefficients a k are lin- 
ear combinations of the system parameters a and 3 and are of 
the form: 

a k = b k a + c k^ + d k (2-13) 



where b k , c k , and d k are real constants. 
Substitution of equation (2-13) into equation 

B^a + c^B + = 0 

B 2 a + C 2 B + D 2 = 0 



(2-12) yields: 



(2-14) 



where 



B n = 



n 

I 

k=0 



B_ = 



n 

l 

k=0 



C, = 



n 

I 

k=0 



C. = 



n 

I 

k=0 



( 1) b k u n U k- 1 ( C) 
( - 1)k Vn U k ( ^ 

(-l) k c. w k U, n (?) 
k n k-1 * 

(-l) k c, a) k U, (?) 
k n k 



(2-15) 
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(2-15) 



n 



D 1 “ I 

1 k=0 



( - 1)k Vn°k-l ( 5> 



n 




(-l) k d, to k U, U) 
k n k 



Solution of equation (2-14) by Cramer's Rule yields the 
following parameter plane solution equations : 



a 



C 1 D 2~ C 2 D 1 

B 1 C 2 -B 2 C 1 
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B 2 D 1 B 1 D 2 
B 1 C 2- B 2 C 1 



(2-16) 



The functional dependence of , B 2 , , C 2 / D^, D 2 on 

£ and w will be omitted for simplicity in this report. A 
list of the algebraic form of the Chebyshev functions of the 
first kind U^(^) is given in Appendix I. Appendix II gives 
a table of U^(?) values for selected values of £ . 

The parameter plane, as defined by Siljak (2) is a rec- 
tangular coordinate plot with a as the abscissa and 6 as the 
ordinate. Equation (2-14) permits the mapping of points, 
excepting real axis points, from the s-plane into the param- 
eter plane. 

For example, C may be fixed <;=<;, and varying to from 

s n 

zero to infinity results in a constant zeta curve in the 
parameter plane. This curve specifies the a, 6 pairs which 
will cause equation ( 2 - 1 ) to have a pair of complex roots 
with damping ratio C = C . Similarly fixing <0 = u> ns and 

varying 5 from zero to one will produce a curve on the param 
eter plane which will specify the a, 8 pairs for which a pair 
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of complex roots in equation (2-1) have a constant natural 
frequency, co . It should be noted that this mapping is not 
conformal since a pair of complex roots on the s-plane maps 
into one point on the parameter plane and, as will be seen, 
a real axis point on the s-plane maps into a straight line 
on the parameter plane. 

For real axis points in the s-plane note that C = 1 and 
letting to n = -a, the equation for s becomes s = -a. The 
characteristic equation becomes: 

n , 

F (s) = l a. (-a) K = 0 . (2-17) 

k=0 K 



Substitution of equation (2-13) into equation (2-17) yields 



n 



• F (s ) = l (b k a+c k 6+d k ) (-a) = 0 



k=0 



Simplifying equation (2-18) yields: 

aB (a) + 0C(a) + D(a) = 0 

where : 

n 

I 

k=0 



B (a) = l (-l) k b k a k 



(2-18) 



(2-19) 



n 



C(a) = l (-1) k c. a k 
k=0 K 



( 2 - 20 ) 



n 



D(ct) = l (-1) k d, a k 
k=0 K 



For a fixed value of a, B(a) , C (a) , and D(a) are constant and 
equation (2-19) in the equation of a straight line on the pa- 
rameter plane. The a,B pairs corresponding to this line cause 
equation (2-1) to have a real root at s = -a. 
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B. MAPPING OF THE PARAMETER PLANE 



The mapping of the s-plane into the parameter plane per- 
mits a designer to choose or adjust parameters so that charac- 
teristic equation roots lie within specified areas of the 
s-plane . 

The rules and graphical techniques for the mapping proc- 
ess are discussed in detail by Siljak (2) and Thaler (13) . 

The major mapping rules and interpretation of the result- 
ing curves are illustrated in the following example. 

Consider a linear fourth order system with the follow- 
ing characteristic equation: 

F(s) = s 4 + (ct+15) s 3 + (15a + 50)s 2 + (50a)s + 6006 = 0 . 

Figure II-l shows the mapping of the characteristic equa- 
tion in the parameter plane. Inspection of Figure II-l re- 
sults in the following observations : 

1) A pair of complex conjugate points on the s- 
plane maps into a single point on the param- 
eter plane. 

2) A real axis point on the s-plane maps into a 
straight line on the parameter plane. 

3) The stable area on the parameter plane, for 
this polynomial is defined by the positive a 
axis, the zeta equal zero curve, and infinity. 

4) Complex roots of the polynomial, for specific 
values of a and 6 are located in the region 
on the parameter plane defined by the zeta 
equal zero and zeta equal one curves. 
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For example, point is at the intersection 
of the u n = 2.0 and s = .4 curves therefore 
the characteristic equation has roots at: 

s = -(.4) (2.0) ± j2.0/l- (.4) 2 
= -0.8 ± jl.83 . 

5) The remaining portion of the stable area 
defined by the zeta equal one curve, the 
positive a axis, and infinity is a real 
root area. For example point is at the 
intersection of four real root lines giving 
roots of: 

s = -1.0, -2.0, -8.0, and -9.0 



15 



006 






004 



x scale 
y scale 



1.0 units/in 

•10 units/in 

4 ^ 

s ( a+15) s + 



(15a+50) s + 50as 



Figure IZ-1 Parameter Piane Example 






6008 = 0 
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III. SINGULAR LINES ON THE PARAMETER PLANE 



The parameter plane developed by Siljak (2) and the com- 
puter program written by Nutting (7) solved for all of the 
roots of the characteristic equation according to then exist- 
ing theory. However, Bowie (8) showed that for some systems 
the existing parameter plane method did not solve for all of 
the roots of the characteristic equation. Bowie found that 
these undetermined roots were actually located on a constant 
zeta-constant natural frequency line which has been defined 
as a singular line. 

A. MATHEMATICAL BASIS OF SINGULAR LINES 

As pointed out in Section II the solution of the parameter 
plane values (equation 2-16) was carried out by the applica- 
tion of Cramer's Rule to equation (2-14). The application of 
Cramer's Rule to equation (2-14) is valid only when the aug- 
mented coefficient matrix of this equation is non singular or 
of rank equal two. Bowie showed that when the augmented 
coefficient matrix of equation (2-14) is singular and the 
rank is equal to one, a condition exists such that 8 becomes 
a linear function of a and a singular line results. In or- 
der for the augmented coefficient matrix to be singular and 
of rank equal one the following conditions must prevail: 
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0 



B 1 C 1 
B 2 C 2 



= B 1 C 2 " B 2 C 1 = 



(3-1) 



- D 1 C 1 
-D 2 C 2 



- C 1 D 2 C 2 D i “ 0 



(3-2) 



B 1 -°1 



B 2 -D 2 



= - Bi d 2 = 0 



(3-3) 



For example, let the solution of equations (2-15) give 
the following values for a specific £ , w n pair: 

B 1 = 2 B 2 = 4 

C 1 = 5 C 2 = 10 

D 1 = 6 D 2 = 12 

Equations (3-1, 3-3) are satisfied and equations (2-14) 
become : 



2a + 56 + 6=0 

4a + 106 + 12 = 0 . (3-4) 

The resulting relationship between a and 6 is: 

D 2a+6 



which is a straight line on the parameter plane. This line 
is a locus of a, 6 pairs which will cause the system charac- 
teristic equation to have a pair of fixed complex roots with 
values of t, and u> n corresponding to those used in solving 
equation (2-15) above. 
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B. METHODS FOR SHOWING EXISTENCE OF SINGULAR LINES 



In general one or more singular lines may be shown to ex- 
ist for a particular system by the application of equations 
(3-1, 3-3) however, this approach requires the substitution 
of specific £ , to n pairs into equation (2-15). Since there 
are an infinite number of Z , , to n pairs which might satisfy 
equations (3-1, 3-3), it is apparent that additional criteria 
for the possible £ , w pairs must be established. 



Bowie derived the following method for finding singular 
lines on the parameter plane (8). Proof of the following 
equations will not be shown but may be found in reference (8). 

The expansion of B^C 2 - B 2 C 1 s ^ own as equation (3-5) 
on the following page, or in compact notation 



Method I 



n 



B..C- - B 0 C. = 7 ce . . = 0 

12 2 1 . L . i] 

i=l J 

j-1 



(3-6) 



where : 



CE = [ce. .] is an (nxn) 
ID 



upper trangular array whose elements are: 




^ . ( C) ^ '*’ + ^ ^ for all i+j < n+2 

3 n J 



0 



for all i+j > n+2 



(3-7) 



Similarly : 



n 



C n D 0 - C-Dt = J ce. . = 0 

12 21 . L . i] 

1=1 J 



(3-8) 



j = l 
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( 3 - 5 ) 



where : 



CE = [ce. .] is an (nxn) 
upper triangular array whose elements are: 



[c. , d . . 1 - c . , . ,d. J(-l) 1 1 U. (C)aj 2l+j 2 for all i+j < n+2 

l-l i+j-1 i+j-1 l-l 1 n J 



and 



n 



B 1 D 2 “ B 2 



D. = - Y ce. . = 0 
1 i£l ^ 

j=l 



where 



for all i+j >_ n+2 
(3-9) 

(3-10) 



CE = [ce^j] is an (nxn) 
upper triangular array whose elements are 



[b. . d . . . -b . , . , d . . ] (-1) 1_1 U. U)co 2l+j 2 

l-l 1+3-I l+i-l l-l in 



i+D 



0 



for all i+j < n+2 
for all i+j > n+2 



(3-11) 



Equations (3-6), (3-8) and (3-10) give three polynomials 

in to and t, of order 2N-1 ' in to and N-l in For simple sys- 

n / n 

terns it is possible to solve for unique values of £ and co n 
which satisfy equations (3-6), (3-8) and (3-10). For more 

complex systems it is necessary to set C = ? s in one of the 
equations, solve for to n , and check the resulting values in 
the two remaining equations. This approach is best handled 
in a digital computer. 
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Example I 

Let the system characteristic equation be: 



3 2 

F(s) = s + as + Bs + 2.5B - 12.5 = 0 . 

Solving equations (3-6), (3-8), and (3-10) the follow- 

ing is obtained: 

B 1 C 2 - B 2 C 1 * 5 '“n - “n - 0 

5 2 

B 0 D , - B n D 0 = -co + 2 5 Coo = 0 
2 1 12 n n 

C,D 0 - C-D. = -2 w 4 ? + 2.5 o) 3 (4c 2 -1) + 12. 5u = 0 . 

1 z z ± n n n 

Solving the above for ? and w the following unique val- 
ues are obtained: 






1 _ 

/5 



co 



n 



/5 



Applying these values to equation (2-14) results in: 
5a - 2 . 5B = -2.5 



and 



2/5 a - /5 B = -/5 
and the singular line relation is: 

B = 2a + 1 



Method II 

The following approach to the singular line problem was 
developed by Cadena (10) . 

Assuming that a system has a singular line for some 
a> n pair, there must be a pair of complex roots which are 
fixed for all a,B pairs in the parameter plane. Therefore, 
let a pair of complex roots of the characteristic equation 
be given by: 



22 



(3-12) 



s = -0 ± j/o ) 2 -a 2 
J n 

where a = After multiplication and manipulation the 

equation (3.12) becomes: 

s 2 + 2a + w 2 = 0 . (3-13) 

n 

If this pair of roots is divided into the characteristic 
equation there results a quotient which is in effect a re- 
duced characteristic equation which has as roots the remain- 
ing roots of the original characteristic equation. Also 
there is, in general, a remainder in the form: 

1 k 

£ s [f (a, 3 ,a ,w ) ] . (3-14) 

k=0 n 

If the complex roots, chosen as the divisor, are to be 
roots of the original characteristic equation, the coeffi- 
cients of equation (3-14) must be independently equal to 
zero. Setting t;hese coefficients equal to zero with the re- 
striction that 13 is a linear function of a results in a 
system of polynomials which is the same as the system devel- 
oped in Method I above. 

Example 1 . 1 

Dividing the characteristic equation: 

F (s) = s 3 + as 2 + 3s + 2.53 - 12.5 = 0 

(from Example I above) by: 

2 2 
s + 2as + w 

n 

results in a quotient or reduced characteristic equation of: 

s + a - 2 = 0 
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and a remainder of 

(6-o) 2 -2aa+4a 2 ) s + (2 . 5B-aio 2 +2aco 2 -12 . 5) = 0 . 
n n n 

Since both terms of the remainder must be equal to zero 
and since 6 is a linear function of a the derivative of B with 
respect to a must be the same for both terms. From the first 
term: 



dB 

da 



2a 



and from the second term 



dB 

da 



2w 



n 



Equating the derivatives : 

2 c 
a) = 5a 
n 

Substitution of a = £io yields : 

n 1 

“n = 55 • 

Returning to the first remainder term, solving for B and sub- 
stituting for B in the second remainder term results in: 

a (10a-2oj 2 ) + 5 to 2 - 20a 2 - 25 + 4aco 2 = 0 . 

n n n 

Since the above equation must be satisfied for all a it must 
be true that the a coefficient is equal to zero and that the 
remaining constant terms sum to zero. From the a coefficient: 



2 c 
lo = 5a 
n 

which was derived above using derivatives. 
Summing the constant terms : 



5 to 2 - 20a 2 - 25 + 4aoo 2 = 
n n 
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Returning to the first remainder term, solving for a in 

terms of 8, substituting into the second remainder term and 

setting the sum of the constant terms equal to zero results in 

4 OC 
co = 25a 
n 

Summarizing the above relationships and letting a = CU) n 
results in: 



to = 5c 
n 

w n = Tc t5^(4? 2 -l)+25] 



co = 2 5C 
n 



Comparing these results with the results of Example I 
above shows that the singular line relationships are identical 
Solving the polynomials for C and co results in: 

c -i- 

/5 

co = /5 
n 

a = 1 

Replacing these values in one remainder term provides 
the required a, 6 relationship for the singular line: 

8 = 2a + 1 



The third root is obtained by substitution of o into the 
quotient term giving: 

s = -a + 2 
Method III 

Karmarker derived the following method for showing the 
existence of singular lines in the parameter plane (11) . 
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The system characteristic equation may be written as: 



n , n-1 

s + a s 
n 



+ a n s + a =0 
1 o 



(3-15) 



where : 



a k ■ V + c k B + d k 



(3-16) 

Let one pair of complex roots for equation (3-15) be: 



s = -a ± j/t o 2 -a 2 
n 



(3-17) 

Now equation (3-15) may be factored into two polynomials: 

(3-18) 

which is obtained from equation (3-17) and: 



s^ + 2as + = 0 

n 



n-2 n-3 

s + s 



x 0 s + x „ = 0 
n-3 n-2 



(3-19) 



which is the reduced characteristic equation. 

Multiplying equation (3-18 and (3-19) equating to the 
coefficient in equation (3-15) and putting the matrix nota- 
tion results in: 



where 



A = 





[A] [B] 


= [C] 






n-1 


C T 
n-1 


-1 


0 


0 ... 


0 


n-2 


0 
3 

1 

to 


-2a 


-1 


0 


• 


• 


• 


2 

-03 

n 


-2a 


-1 


• 


• 


• 


0 


2 

-CO 

n 


-2a 


• 


• 


• 


• 


0 


2 

-w 

n 


• 



0 

0 



0 

0 



0 

0 

0 



-1 

-2a 

2 

-a) 

n 



(3-21) 
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; a 




20 - d n-l 


6 




to 2 - d - 
n n-2 


X 1 




-d _ 
n- 3 


x 2 




• 


• 


c = 


• 


• 




• 


• 

X ~ 

n- 3 




• 

- d i 


x 0 
n-2 




- d o 



For singular lines to exist 6 must be a linear function 
of a. This criterion plus the techniques of linear algebra 
result in the following singular line existence theorem: 
Singular Line Existence Theorem 

A. All determinants of order n that can be formed 
from the augmented matrix [A,C] are zero; (n+1) such 
determinants may be formed. Note that n of these 
form an independent set. 

B. At least one of the determinants of order (n-1), 
that can be formed from the augmented matrix [A,C] 
must be non-zero. There are 

n+1 n 



n-1 n-1 



such determinants. 
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The conditions that must exist for singular lines may be 



found using the existence theorem as shown in the following 
example . 

Example III 

Consider the characteristic equation used in Example II: 
F ( s ) = s 3 + as 2 + 6s + 2.56 - 12.5 = 0 





10-1 




a 




2a 


[A] = 


0 1 -2a 


[B] = 


6 


[C] = 


2 

to 

n 




0 2.5 -to 2 

n 




X 




12.5 



The (n+1) determinants of order n which can be formed 
from the augmented matrix [A,C] are: 

















1 


0 


-1 




1 


0 


2a 


0 


1 


-2a 




0 


1 


2 

to 

n 


0 


2.5 


2 

"“n 




0 


2.5 


12.5 



1 2a -1 




2a 0 -1 


0 

CM 

1 

CM 

3 

o 




to 2 1 -2a 


n 




n 


0 2.5 -to 2 




12.5 2.5 -to 2 


n 




n 



Solving the above determinants and equating each to zero 

results in the following: 

2 c 
to = 5a 
n 

to^ = 25a 
n 
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2 _ 10a 2 +12 . 5 
w n 2a+2 . 5 




Solving for a and to : 

a = 1 

co = /5 
n 

It should be noted that replacing a with Cco n in the first 
three equations given above results in the same set of equa- 
tions as given in Examples I and II. The last equation is re- 
dundant. 

Substituting for a and to in equation (3-20) and solving 
results in: 

6 = 2a + 1 
r^ = x = a - 2 

Again, this is the same result obtained in Examples I 
and II above. 



Method IV 

The following parameter plane equation and singular line 
theory was derived by Randolph (12). 

A second derivation of the parameter plane transformation 
which is somewhat different than that described in Section II 
is described below. 

Let the system characteristic equation be written as: 

n V 

F(s) = ) a,s K = 0 (3-22) 

k=0 k 
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where 



(3-23) 



a k ■ V + c k 6 + d k • 

The complex variable s may be written as: 

s k = to k e^ k ® = to k [cos k9 + jsin k9] (3-24) 

where 9 is the angle from the positive real axis to a point 
in the s-plane and to n is the natural frequency of the system. 

Substituting equation (3-24 and (3-23) into equation 
(3-22) results in: 

n k 

F(s) = 7 [b.ot+c.B+d, ]to [cos k9+jsin k9] = 0 . (3-25) 

k=0 K K K n 



Letting the real and imaginary terms go independently to 
zero results in 



B + C 1 6 + D 1 = 0 



B2 a + C2B + D2 = 0 



(3-26) 



where 



n , 

B.. = 1 b, to cos k9 

1 k=0 X n 



n 



B 0 = y b, to sin k9 
2 k£o k n 



V k 

C. = V c, to cos k9 
1 k£o k n 



n 



C„ = y c, to k sin k9 
2 k=0 k n 



n 



D 1 = l <3 to cos k9 
k=0 K n 



n 



D 0 = J d, to k sin k9 . (3-27) 

Z k=0 K n 



It should be noted that the coefficients of equation 
(3-26) are not the same as described in Section II above. 

From the theory of linear algebra it may be shown that 
two functions, f 1 (a,B) and f 2 (a,S) are linearly dependent if 
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0 for all 



and only if K^f^(a,B) + K^f^(a,3) + I< 2 f 2 (oi,B) = 
a, 6 in a domain common to f^ and f 2 for some and K 2 not 
both zero. 

Let: 

f^(a,3) = B^a + C^3 + 

f 2 (a, 6 ) = B 2 a + C 2 3 + D 2 . (3-28) 

Then a necessary and sufficient condition that f^(a,B) 
and f 2 (ct, 6 ) be linearly dependent is that there exists a con- 
stant M such that: 



b 2 = m B;l 



C 2 = MC 1 



d 2 = m C;l 



(3-29) 



The conditions that must be satisfied in order for a 
singular line to exist, may be obtained as follows: 

(a) Write equation (3-27) . 

(b) Relate these equations using equations (3-29) . 

(c) Solve for M, 0, and . 

(d) The resulting singular line equation is obtained 
from equation (3-26) 

Example IV 

Using the same characteristic equation as in previous 
examples : 

F ( s ) = s 3 + as 2 + 3 s + 2.53 - 12.5 = 0 
Equations (3-29) may be written as 

MB- = M[o ) 2 cos20] = co 2 sin20 = B„ 
in n z 
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MC, = M[2.5+to cosB] = to sin0 = C 0 

1 n n Z 

MD. = M [-12 . 5+to 3 cos30] = to 3 sin30 = D 

i n n 

Solving for M, 0, and to n results in: 




0 = cos (- — ) 
/5 

to = /5 

n 



Since : 

0 = cos ^(-£) 

This results in a complex pair of roots at: 

s = -1 ± j 2 

Returning to equation (3-26) and solving for the singu- 
lar line equation results in: 

6 = a - 2 

Dividing the characteristic equation by the complex 
roots results in the third root of the system: 

r = a - 2 

These results are the same as those of previous examples. 

The preceding portion of this section developed methods 
for showing the existence of singular lines and the illus- 
trative examples indicated procedures for finding the neces- 
sary conditions for singular lines using the existence the- 
orems. These procedures work well for simple systems, but 
become complicated for the higher order systems since, in 
general, the order of ton is 2 (n-1) and the order of a is 
n-1 in the resulting equations. 
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The following procedure provides a simple technique for 
finding singular line characteristic equations. Manipulation 
of the resulting characteristic equations will in general pro- 
vide a singular line plant. 

Let a system characteristic equation be written as: 

n , 

F (s ) = l s a, = 0 (3-30) 

k=0 K 



where : 

a k = ab k + $ c k + d k (3-31) 

Assume that F(s) is a singular line characteristic 

equation for some £=£,<») = w . Then equation (3-30) 

s n ns 

may be rewritten in factored form as: 

F (s) = (s^+2t u) +to^ ) [ (s+r, ) (s+r 9 ) .... (s+r 0 ) ] = 0 
s ns ns 1 z n-z 

(3-32) 

In general the roots r^ , , .... ^ n _ 2 ' are functions 

of a and 3. However for the singular line case a and 8 
must appear in linear terms such that the order of a and 8 
is one and a8 products do not appear. Therefore only one 
root of r-^^ • • ♦ • r n _ 2 ' ma Y contain terms in a and 8 and this 
root must be real. The remaining roots may be fixed at any 
desired value. If the resulting mobile real root contains 
only a as a variable a^ will be a function of a. A singular 
line may now be obtained by substituting 8 as a linear func- 
tion of a into one or more of the resulting a^ terms. 

For example suppose that for a third order system a 
singular line is desired such that the roots are: 
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s = -1 ± j2 



s = a - 2 

Substituting into equation (3-32) yields: 

F(s) = (s 2 +2s+5) (s+a-2) = s 3 + as 2 + (2a+l)s + 5a - 10 = 0 . 

(3-33) 

Let 3 = 2a + 1 and substitution into equation (3-33) results 
in : 

F ( s) = s 3 + as 2 + 6s + 2.56 - 12.5 = 0 . (3-34) 

Equation (3-34) is a singular line characteristic equa- 
tion for 6 = 2a + 1 as shown in previous examples. 

The procedure indicated above results in a unique sin- 
gular line for a given characteristic equation. A more gen- 
eral case is examined by Bowie (8) in which there are, 
apparently, an infinite number of singular lines for a given 
characteristic equation. The system examined is a complex 
cross-coupled type with powers of s to the sixth power occur- 
ring in the characteristic equation. There does not appear 
to be a tractable method for producing similar systems at 
the present time. 
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IV. EXAMPLES OF SYSTEMS WITH SINGULAR LINES 



This section contains block diagrams, characteristic 
curves, and step response characteristics for systems which 
have singular lines. The systems are grouped by the nature 
of their characteristic equations. 

A. SYSTEMS WITH ONE SINGULAR LINE AND ONE MOBILE ROOT 



Ex. 


Characteristic 

Equation 


Roots 


Singular 

Line 


Figs . 


1 


3 2 

s +as +(2a+l)s+6=0 


s=-l± j 2 
s=— (a-2) 


6=5 (a-2 ) 


IV-1 

IV-2 


2 


s 3 +as 3 +6s+2 .56-100=0 


s=-l± j 2 
s=- (a-2) 


6=2a+l 


IV- 3 
IV- 4 


3 


s 3 +as 3 + (46+5) s+ 6=0 


s=-l± j 2 
s=- (a-2) 


6=5 (a-2) 


IV- 5 
IV- 6 


4 


s 4 + (a+10) s 3 + (12a+l) s 3 
+ 25as+10B=0 


s=-l± j 2 
s=- (a-2) 
s=-10 


6=5 (a-2) 


IV- 7 
IV- 8 


5 


s 4 + (a+10) s 3 + (12a+l) s 3 
+ 6s+2B-100=0 


s=-l± j 2 
s=- (a-2) 
s=-10 


B=25a 


IV- 9 
IV-10 


6 


s 4 + (a+10 ) s 3 + ( 12a+l) s 2 
+ (56+50) s+106=0 


s=-l± j 2 
s=- (a-2) 
s=-10 


6=5 (a-2) 


IV- 11 
IV-12 



As a is increased the mobile root moves further from the 
origin of the s-plane and the effect of this root becomes 
negligible for values of a greater than 20. This trend is 
shown by the step response characteristics for the respective 
systems . 
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The movement of the mobile root along the real axis of 
the s -plane maps along the singular line in the parameter 
plane; this is shown by substituting x = -(a-2) into the sys- 
tem characteristic equation. In the fourth order systems, 
the fixed real root, s = -10, also maps along the singular 
line . 
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Figure IV-1 Step Response Curves - Example 




5T C 2T0 OTO 




x scale = 20 units/in 
y scale = 20 units/in 

s^+as^+(2a+l)s+6=0 

Figure IV-2 Parameter Plane Curves - Example 1 
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Figure IV-3 Step Response Curves - Example 




x scale = 20 units/in 
y scale = 20 units/in 

s 3 +as 2 +6s+2. 56-100=0 

Figure IV-4 Parameter Plane Curves - Example 2 
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Figure IV-5 Step Response Curves - Example 



6T0 2T0 0 TO 800 900 600 200 000 




x scale = 20 units/in 
y scale = 20 units/in 

s 3 +as 2 + (46+5) s+$=0 

Figure IV-6 Parameter Plane Curves - Example 3 
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Figure IV- 7 Step Response Curves - Example 



900 500 600 V30 200 TOO OOO 




x scale = 3 units/in 

y scale = 10 units/in 

s 4 +(a+10) s 3 + (12a+l) s 2 +25as+106=0 

Figure IV-8 Parameter Plane Curves - Example 4 
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STO 2 TO 600 900 EDO 000 
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Figure IV-9 Step Response Curves - Example 





x scale = 3 units/in 

y scale = 10 units/in 

s 4 + (ct+10) s 3 + (12a+l) s 2 +Bs+2B-100=0 
Figure IV-10 Parameter Plane Curves - Example 5 
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S20 020 STO OTO SOO 000 




x scale = 2 units/in 
y scale = 5 units/in 

s 4 + (a+10) s 3 +(12a+l) s 2 + (56+50) s+106=0 
Figure IV-12 Parameter Plane Curves - Example 6 
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SYSTEMS WITH POLE - ZERO CANCELLATION 
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Figure IV-13 Step Response Curves - Example 



6 TO 2T0 OTO 800 900 WO 200 000 




x scale = 20 units/in 
y scale = 20 units/in 

s 3 + (8+35) s 2 +( 308+550) s+200 (a+6+5) =0 

Figure IV-14 Parameter Plane Curves - Example 7 
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Figure IV-15 Srep Response Curves - Example 




x scale = 20 units/in 
y scale = 20 units/in 

s 4 +(B + 40) s 3 +40 (6+10) s 2 + 400 (a+7. 5) s+3000ct=0 
Figure IV-16 Parameter Plane Curves - Example 8 
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C. SYSTEM WITH LOCKED ROOTS 



The system shown in Figure IV-17 has the following char- 
acteristic equation: 

s 3 + (a+453+22) s 2 + (22a+9903+45) s + 100 = 0 . 

Setting a = -453 forces the system to have the following 
roots : 

s = -20 
s = -1 ± j 2 

There are no parameter plane curves for this system since 
the parameter consists of an infinite number of singular 
lines . 
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X >1 
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Figure IV-17 Step Response Curves - Example 



V. SELF-ADAPTIVE CONTROL USING SINGULAR LINES 



One important application of singular line systems is 
that of self-adaptive control. For example, a parameter 
which is a function of environmental conditions could be 
compensated for by linearly varying a second parameter. 

This situation might occur in an aircraft where the dynamic 
flight characteristics are a function of altitude. 

Consider Figure V-l, below, which shows a singular line 
for Z = C s and w n = u) ns on the parameter plane. 3 is the 
varying plane parameter and a is an adjustable parameter in 
the controller. The initial operating condition is on the 
singular line at point so that dominant complex roots 
are at 



s = ~Z w ± jw / 1~Z 2 
s ns J ns s 



When 3 is varied by an amount A3 the operating point is 
changed to point I >2 by changing a by an amount, Aa. The sys- 
tem maintains roots at 



s = -£ a) ± jui / 1-Z 2 
s ns J ns s 

The other characteristic roots may be changed, but for rea- 
sonable changes in a and 3 root dominance is retained. 
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The system in Figure V-2 is an example of a singular line 
plant (pole-zero cancellation) which has been designed for 
self-adaptive control by a procedure developed by R. A. 
Desrosiers (13). Explanation of the various control func- 
tions will not be explained but may be found in reference 
(13) . Figure V-2 also shows the step response of the system 
with a 12% step change of 6 at t - .09 sec. As indicated, 
the adapted plant and model responses coincide. (The model 
is the basic plant with no parameter changes.) 
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VI. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER WORK 



Feedback control systems containing singular lines , while 
occurring in a limited number of cases, provide a tool for 
the design engineer working in self-adaptive control and 
other areas. Methods for predicting the existence of singu- 
lar lines have been shown with examples given in the use of 
these methods. As indicated in previous sections of this 
report, the equations resulting from the application of known 
existence techniques are complex and do not, in general, pro- 
vide a tractable means of finding singular line systems. 

From the standpoint of the design engineer, it would be 
desirable to formulate a general set of rules which would 
provide a specific system with singular lines. Research to 
date does not indicate the existence of such a general set 
of rules. Further development is needed in this area. 

A second approach might be to design a number of singu- 
lar line systems using characteristic equations as starting 
points. A set of typical singular line systems would then 
be made available to the design engineer for correlation be- 
tween a specific system and possible singular line charac- 
teristics . 

Another area for future work might be to develop further 
theorems for singular lines on the parameter plane. For ex- 
ample, all real root lines for singular line systems were 
found to lie on the singular line. Further study may prove 
that this is true in general. 
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APPENDIX I 



TABLE OF THE CHEBYSHEV FUNCTIONS OF THE SECOND KIND 

The Chebyshev functions of the second kind, are de- 

fined by the recursion relation: 

U k+l U) " 2 ^ U k ( ^ + U k-l U) = 0 

where : 





V ? > 


= 0 




U 1 (?) 


= 1 . 


functions 


for k = -1 , . 


. . , 10 are 


k 


Vo 




-1 


-1 




0 


0 




1 


1 




2 


2 £ 




3 


4c 2 


- 1 


4 


8 , 3 


- 4C 


5 


16? 4 


- 12C 2 + 1 


6 


32 C 5 


- 32C 3 + 6 C 


7 


64? 6 


- 80c 4 + 24c 2 - 1 


8 


128? 7 


5 3 

- 192c + 80C “ 8c 


9 


256C 8 


- 448C 6 + 240C 4 - 40C 2 + 1 


10 


9 

512C 


- 1024C 7 + 672C 5 - 160C 3 + 10c 
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